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Upper limit on the gas density in the f3 Pictoris system 
On the effect of gas drag on the dust dynamics 
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Abstract. We investigate in this paper the effect of gas drag on the dynamics of the dust particles in the edge-on 
(5 Pictoris disc in order to derive an upper limit on the mass of gas in this system. Our study is motivated by the 
large uncertainties on the amount of gas in the j3 Pictoris disc currently found in the literature. The dust particles 
are assumed to originate from a colliding annulus of planetesimals peaked around 100 AU from the central star as 
proposed by Augereau et al. (2001). We consider the various gas densities that have been inferred from independent 
observing techniques and we discuss their impact on the dust dynamics and on the disc profile in scattered light 
along the midplane. We show that the observed scattered light profile of the disc cannot be properly reproduced 
if the hygrogen gas number density at 117 AU exceeds 10 4 cm -3 . This corresponds to an upper limit on the total 
gas mass of about 0.4 assuming the gas density profile inferred by Brandeker et al. (2004) and thus to a gas to 
dust mass ratio smaller than 1. Our approach therefore provides an independent diagnostic for a gas depletion in 
the f3 Pictoris system relative to the dust disc. Such an approach could also be used to constrain the gas content 
of recently identified systems like the edge-on disc around AUMic. 
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1. Introduction 



planetary systems - stars: j3 Pictoris- planetary systems: formation - planets and satellites: 



The f3 Pictoris system is still one the best known ex- 
amples of a debris disc surrounding a young Main 
Sequence star. The dust component of the almost edge- 
on disc has been thoroughly observationally investigated, 
in thermal emissi on as well as in scattered light (e.g. 
I Ar t vmowiczl 1 1 9 9 7t and references therein). The total dust 
mass is believ ed to range between a f ew lunar masses and 
0.5 Mm (e.g. IZuckerman fe Becklinl Il993t lArtvmowicj , 
ll997tlDent et al.ll200Cj|) . These estimates lead to collisional 
lifetimes of dust particles which are mu ch shorter than the 
estimated age of t he system (8-20 Myr ; IZuckerman et all 
l200ll ISong et all 120031: iDi Folco et ail. l2004|h which im- 
plies that the dust disc cannot be primordial but has to 
be r eplenished by collisions or evaporation of la rger bodies 
(e.g. lArtvmowic2lll997tlLecavelier et all 119961. Scattered 
light midplane surface b rightness profiles displ ay a sharp 
decrease beyond 120 AU l|Kalas fe Jewittlll995|) where the 
luminosity profile follows a steep radial power law r a with 



-5.5 < a < -4.5 l|Heap et alJ.l2000|) . These outer regions 
of the disc that extend over hundreds of AU are believed to 
be made of small grains primarily produced within 150 AU 
and placed on high eccentricity orbits by radiatio n pres- 
sure l|Lecavelier et allll996tlAugereau et al.U200lj) . 

In contrast to our relatively good knowledge of the 
dust distribution, the gaseous compound is still poorly 
constrained. This is in particular the case for hydrogen, 
which should a priori make up most of the gas disc. So 
far, attempts to directly d e tect h ydrogen have failed, with 
the exception of lThi etal 1 fcOOll) ' who have measured with 
ISO H 2 infrared emission features associated to pure rota- 
tional transitions. The mid- infrared observations are con- 
sistent with more than 50 Mm of warm H2 in the Pictoris 
disc. These results are challenged by the non-detection in 
the UV with the FUSE satellite of absorption lines from 
H2 electronic transition l|Lecavelier et alll200l[) . This non- 
detection places an upper limit on the H2 column density 
three order of magnitude smaller than the value derived 
from the infrared observations. 
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The H2 e mission is unlikely to be of interstellar origin 
according to iLecavelier et al I l|l996h . The ISO aperture 
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been large, the infrared emission could originate from re- 
gions not probed by FUSE if the gas disc is very spatially 
extended. Recent mid-infrared observations with Spitzer 
nevertheless do no confir m the 17.03 5 umH 2 emission fea- 
ture observed with ISO l|Chen et all (2004). The gas mass 
upper limit based on the Spitzer results is 11 M® assuming 
a gas temperature of 110 K. 

Attempts have also been made to indirectly recon- 
struct the hydrogen abundance from other more accu- 
rately observed species. Absorption lines in the stel- 
lar spectrum have indeed reveale d the prese nce of 

numerous met a llic e l ements (e.g. Hobbs et all Il985l 
lLagrange et al 1 Il998t iRoberge et all 1200(1) . From early 
estimates, rescaling to solar abundances lead to col- 
umn densities of NCR) = 10 18 -4 1 20 cm- 2 <|Hobbs et all 
Il985[) . iKamp et alJ l)2003() performed emission line calcu- 
lations for different gas tracers (CO but also C and C + ) in 
the disc. They conclude that the total (3 Pictoris gas mass 
(presumably mostly made of hydrogen) would l ie bet ween 
0.2 and 4M©. More recently. IBrandeker et alJ l(2004h . us- 
ing the VLT/UVES instrument, spatially resolved the gas 
disc and observed numerous emission lines, mainly from 
Fel, Call and Nal which was observed as far as 323 AU 
from the star. Estimating a ionisation structure for the 
disc and assuming solar abundances, they found a column 
density for atomic hydrogen o f JV(HI) = 8 10 18 cm~ 2 , con- 
sistent with the upper limit of lFreudling et ail l|l995|) . and 
iV(H2) = 3 10 18 cm~ 2 for mole cular hydrogen, consistent 
with the lLecavelier et all l|200lj) upper limit. However, the 
inferred total mass of the gas disc, ~ 0.1Af 8 , is incon- 
sistent wi t h the 50 M® gas mass measured with ISO by 
iThi et al.l (|200lh - 

Nevertheless, independent considerations might argue 
in favour of a more massive gas disc. One puzzling re- 
sult deduced from the observed absorption lines is in- 
deed that all ion species seem to be on keplerian orbits , 
with low relative radial velocities l)Ok>fsson et all boOll) . 
This is in strong contradiction with the fact that many 
of these elements should be rapidly ejected on unbound 
orbits by strong stellar radiation pressure (/3 Pictoris is an 
A5V star). Several explanations have been proposed to 
explain this contradiction. One of them is that gaseous 
friction due to an unseen gaseous braking agen t, possi- 
bly h ydrogen, should damp the outw ard outflow (iL jsgjuiJ, 
2003). With calculations based on the lLiseau et alJ|l999|) 
model, IBrandeker et alJ l|2004l) estimated that a high den- 
sity gas disc of ~ 50 M® could lead to radial velocities 
compatible with the observed ones provided the disc is 
metal depleted. 

There are thus still large uncertainties regarding the 
amount of gas in the (3 Pictoris disc. We propose here to 
address this issue by looking at the effect of gaseous fric- 
tion on the dynamics of the dust grains observed in scat- 
tered light. This approach is used to constrain the gas den- 
sity in the disc and its r adial distribution . We b ase our ap- 
proach on the results of lAugereau et al I ll200lt) . who have 
shown that the observed dust disc could be successfully 
model if the grains are produced by a distribution of par- 



ent bodies within 150 AU from the central star and placed 
on orbits with large semi-major axis (a) and high eccen- 
tricities (e) by radiation pressure. The important point is 
that this satisfying fit of the surface brightness profile was 
obtained for a gas free system. 

Our main goal is here to estimate how gas drag might 
affect these results, and in particular which gas disc densi- 
ties and distributions yield dust density profiles still com- 
patible with observations and which ones don't and might 
thus be ruled out. Our approach is much in the spirit o f 
the remarkable pioneering work of lLecavelier et aD fl996). 
who explored the scattered light profile for different par- 
ent bodies prescriptions and even included gas drag, but 
considered an academic gas disc prescription (p g cx 1/r) 
and restricted their approach to (3 < 1 particles (note 
that the behavi o ur of j3 > 1 particles was addressed by 
iLecavelier et alJ l|l998|) in the case of BD+31643, a very 
different stellar environment). This neglection of the /? > 1 
particles might be too restrictive in the presence of effi- 
cient gas drag which prevents high [3 particles from leaving 
the system as quickly as in a gas free disc (see section 2.2 
and Fig-EJ ■ Furthermore, at that time the radial power 
law index for the outer midplane luminosity profile was 
believed to be in the [—3.6, —4.3] range rather than in the 
[—4.5, —5.5] range. 

2. Model 

2.1. Equations of motion 

Our integrator is a classic 4 th order Runge-Kutta which 
test runs (performing mutual comparisons with higher or- 
der algorithms) have proved to be of a satisfying accuracy 
for the present study. The integrator takes into account 
the star's central potential, the radiation pressure force 
and gaseous friction. The radiation pressure on a grain of 
mass m at a distance r from the central star is expressed as 
a function of the gravitational potential trough the usual 
(3 parameter 



- rad 



— f3Fg rav p 



GMm 



(1) 



assuming an optically thin disc in all directions. Here G is 
the gravitational constant and M the mass of the central 
star. We will consider grains larger than 0.1 /j,m radius, 
for which there is an inverse one t o one relation between 
(3 and grain size s l|Artvmowiczl Il988j) and we assume 
(3 ~ 0.5 x (5/um/s) with s expressed in ura. For th e gas 
drag forc e we fo l low iTakeuchi &: Artvmowiczl l|200lj) and 
adopt the iKwokl l|l975|) formalism 



,s 2 (v 2 T + Av 2 ) 1/2 Av 



gas 



(2) 



where p g is the gas density, vt is the gas thermal velocity 
and Av is the difference between the grain velocity v and 
the gas velocity v g expressed as a function of the keplerian 
velocity i>k 



v g = i>k(1 - vf 



(3) 
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Distance from the star in AT." 

Fig. 1. Surface density of the parent bodies disc responsi- 
ble for the reso l ved sc attered light images according to 
lAueereau et alJ l|200lL solid line). The dashed line to- 
gether with the right y axis give the power-law index of 
the profile. 

The parameter v depends on the gas pressure gradient 
dPg/dr and reads 

1 dP g 
rfl 2 p g dr 

2.2. Dust disc 

We assume that the dust particles are produced by par- 
ent bodi es on circular orbits fol lowing the distribution de- 
rived bv lAueereau et all l)200ll) and displayed in Figure^ 
This parent body distribution has been shown to pro- 
duce a spatially extended dust disc that correctly matches 
the resolved scattered light images as well as the long- 
wavelength photometric data. The parent body distribu- 
tion also results in a simulated disc radial brightness pro- 
file at 12 q m consistent with th e N-band thermal images 

I I i I 

obtained by Pant in et alJ 1119971) b evond ~30 AU. But im- 
portantly, lAugereau et alJ Ij2f)f)ll) couldn't reproduce the 
absolute flux that they attribute to an additional popula- 
tion of very small grains with (3 ratios likely larger than 
0.5. Recent mid- i nfrared observations supp ort this picture 
( Okamoto et"afl 120041: iTelesco et aD. 120051) . The adopted 
distribution of parent bodies displayed in Fig.^ peaks in 
the 80-120 AU region and cuts beyond 150 AU (for the 
validity of this assumption in the context of the present 
paper see discussion in section 4). 

The differential size distribution of the produced dust 
is assumed to follow the usua l equilibrium power law dA*" = 
C s~ 3 - 5 ds llnohnanvil . ITflfiflh . Wc then consider an arbi- 
trary number N num of particles satisfying this distribu- 
tion between s m i n = 0.1 /zm ((3 = 25) and s max = 500 /im 
((3 = 0.005). We consider objects up to high (3 because gas 
drag might be very efficient in slowing down the escape of 
these small grains, hence increasing their contribution to 
the luminosity profile (for a more quantitative analysis, 
see section 3). In order to get statistical significance for 



all size ranges despite the steep power law of the size dis- 
tribution, three different runs, each with N nunl = 20000, 
are performed for the 0.005 < (3 < 0.05, 0.05 < (3 < 0.5 
and 0.5 < [3 < 25 populations, which are then recombined 
with the appropriate weights. 

At each moment every particle is assigned a collision 
destruction probability, depending on its size, velocity and 
location in the disc, and has thus a limited lifetime. Due to 
obvious computational constraints, a realistic modeling of 
collisions is here excluded. We use an empirical approach 
where each particle is assigned, at each timestep, a col- 
lisional destruction probability depending on its (3 value, 
its distance r from the central star and its radial velocity 
v r 

/coil = t^- (5) 

tcoll 

where At is the simulation timestep and i co n the collision 
time scale. We approximate t co n by 

( /3q\ a ( n lmAU \ ( r \0-5 y 
tc ° U= {j) {-^-){l00Mj) — tcM " (6) 

where t co \\ ~ 10 4 years is a referenc e collision timescale 
that we take equal to the estimate of lArtvmowiczl l|l997fl 
for a (3q = 0.3 particle at 100 AU from the central star. v ro 
is the radial velocity of the (3q particles at 100 AU and nji 
is the volumic number dust density at distance r. The (3 
(i.e. grain size) dependency of i co ii relies on complex phys- 
ical parameters. How ever, numeric al studies (e.g. F i g. 7 o f 
IWvatt & Dentl l|2002l) and Fig. 12 o flThebault et alll2003|) 
indicate that the coefficient a ranges between 0.3 and 0.5 
for small particles. We take here a — 0.4. The absolute 
value of the collision rates might be rescaled by changing 
the t co \\ parameter. 

For n r in Eq.Elwe as sume t he dust density profile de- 
rived bv lAueereau et aD l)200l|) (see their Fig. 2) from the 
parent bodies profile displayed in Fig.^ A more coher- 
ent approach would have consisted in taking into account 
n r matching the final dust distribution for each individ- 
ual simulation we performed. But this would have added 
a lot of complexity to the problem by introducing a cir- 
cular argument, the final distribution being in principle 
dependent on the assumed collision timescale. We never- 
theless empirically checked the validity of our prescription 
by performing one test run with n r matching the final dust 
distribution obtained in the reference case of Sec. 3.1, and 
ended up with a final luminosity profile which never de- 
parts by more than 10% from the nominal one. Possible er- 
rors due to the inevitable departure between the assumed 
n r prescription and the final dust distribution are thus 
expected to be negligible compared to the impact of the 
parameters explored throughout the paper on the calcu- 
lated brightness profiles. 

Such an empirical prescription is of course but a crude 
approximation of the real collisional behaviour of dust par- 
ticles. To the first order, it provides nevertheless a good 



4 



Thebault and Augereau: Gas density in j3 Pictoris 



indication on how the collision rate relies on the grain size, 
on the grain velocity and position in the disc 1 . 

2.3. Gas disc characteristics 

Several test gas discs, with different profiles and spatial 
extents, are explored. Our refer ence case is the profile de- 
rived bv lBrandeker et all (|2004|) from their observations of 
metallic emission lines, in particular Nal, out to 323 AU. 
From the observed Nal profile they deduced the radial 
density profile of hydrogen nuclei to be 



n(H) = n 



5.3 



(7) 



with r = 117 AU. 

We will explore different no values betwe e n the two 
extreme values considered by iBrandeker et alJ l)2004|) : 

1. no = 2.25 10 3 cm -3 , required if solar abundances are 
assumed, which leads to a total H mass of ~ 0.1M®, 

2. a much higher density: no = 10 6 cm~ 3 that implies 
a significant depart ure from solar abundanc es but re- 
quired, according to IBrandeker et all l)2004|) . if hydro- 
gen is to act as the "braking agent" maintaining all 
high elements on keplerian orbits. 

We will also consider the case where the gas disc is trun- 
cated at 150 AU. Finally, an acade mic case where the gas 
distribution matches the classical lHavashil l)l98lh radial 
profile proportional to r -2 75 is also tested. All initial con- 
ditions are summarized in Tabled 

2.4. Scattered light profile 

The aim of the paper is to use the shape of the scattered 
light surface brightness profile along the disc midplane to 
set an upper limit on the gas density in the Pictoris 
system. We detail in this section the procedure we used to 
calculate the brightness profile. 

All particles are produced at an initial moment to af- 
ter which we let the system dynamically evolve. At regu- 
lar time intervals we compute the brightness profile along 
the disc midplane assuming a perfectly edge-on disc. At 
each projected distance from the star and because the 
disc is assumed to be optically thin, the surface bright- 
ness is obtained by integrating along the line of sight the 
flux scattered by the dust particles. The stellar flux re- 
ceived by a grain is assumed to drop with the square of 
the distance from the star (stellar flux dilution in an opti- 
cally thin environment). Scattering is assumed to be gray 
and isotropic. Gray scattering implies that the scattering- 
cross section is proportional to the geometric cross section. 
For the isotropic assum ption, we refer to the discussion in 
I Augereau et alJ l|200lj) who showed that the anisotropic 



a similar, though simpler, cut-off lifetime prescription act- 
ing as a substit ute for collisional dust de ctruction has been 
recently used bv lDeller fc Maddisor] (120051) 



scattering properties of the grains only affect the regions 
below ~ 80 AU as far as the scattered light radial bright- 
ness profile along the disc midplane is concerned. The final 
profile is progressively obtained by adding these instanta- 
neous profiles. This procedure is stopped when no further 
significant evolution of the total profile is observed at pro- 
jected distances smaller than 500 AU, either because parti- 
cles eventually escape the system or are being collisionally 
destroyed. 

Our procedure is implicitly equivalent to assuming an 
arbitrarily constant dust production rate from parent bod- 
ies. The absolute value of the "real" production rate is 
here an unconstrained parameter that does not affect the 
results since we are only concerned with the shape of the 
surface brightness profile and not its absolute intensity. 
In other terms, only the relative size distribution of the 
produced dust matters in our approach and not the exact 
amount of dust produced by collisions. Our results might 
be rescaled to any absolute flux level through the Co coef- 
ficient because the collisional lifetime of dust in the models 
has been scaled to the observed value rather than being 
dependent on the density of material in the disc. 

3. Results 

For sake of clarity, we will focus on three observational 
parameters, namely the two flux ratios in scattered light 



along the disc midplane F A 



-P400/100 — F' 



40/100 



Fa, 



40Au/-FlOOAU 



and 



40oau/-Fiooau and the averaged power law 
index p out of the radial distribution profile measured be- 
tween 150 and 40 AU. The r e ferenc e observational val- 
ues deduced from iHeap et alJ l|2000f) are -F40/100 — 0.3, 



400/100 



40/100 

5 10 -3 (by extrapolating their profile) and 



-5.5 < p ou t < -4.5. 
3.1. No gas 

We first perform an academic run with no gas drag in 
order to provide us with a reference case to which the 
following runs might be compared. 

The synthetic scattered light profile displayed in Fig.O 
is close to the one obtained by lAugereau et al.l ((2001^1 . 



with Ft 



40/100 



0.35 and F £ 



400/100 



5.10 . However, we 



find pout — —4.7 while I Augereau et al.l l)200l|) got p out be- 
tween —5 and —5.5. This small difference has two causes. 
First of all, the presence of high (> 0.5) grains ig nored 
in the dynamical approach of lAugereau et alJ l(200ll) mod- 
erately contributes to the flux in the outer regions (assum- 
ing a distribu tion for the > 0.5 par ticles similar to those 
with ~ 0.5. lAugereau et"afl l)200l[) also found a similar 
trend). Secondly, we c onsider here the collisional lifetime 
of the particles whereas lAugereau et al.l l|200l[) did a phase 
mixing of randomly generated particle orbits. The depen- 
dency of £ co ii with distance to the star tends to increase 
the density of the ~ 0.5 particles which spend a sig- 
nificant amount of time in the outer regions, hence the 
slightly flatter slope beyond 150 AU. However, this —4.7 
slope is still in agreement with the radial power-law index 
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Table 1. Summary of the different gas disc properties assumed in this paper. The densities no are given at a distance 
ro = 117 AU from the star. 







no [cm J ] 


radial profile 


Total Mass [M©] 


comment 


case 


1 


10 b 


Brandcker et al. (2004) 


40 M© 


density required for the "gaseous braking agent" hypothesis 


case 


2 


2.25 10 3 


Brandeker et al. ( 20041 


0.1 M© 


reconstructed from the Nal profile assuming solar abundances 


case 


3 


10 4 


Brandeker et al. ( 20041 


0.4 M© 




case 


4 


10 5 


Brandeker et al. ( 20041 


4M© 




case 


5 


10 6 


Brandeker et al. (20041 


15 M© 


truncated at 150 AU 


case 


6 


7.10 5 


Hayashi (r -2 ' 75 ) 


40 M© 







time (years) 



time (yei 



Fig. 2. Radial position as a function of time, for different values of (3, for one particle initially located at 100 AU: a) 
gas free case b) high gas density case (case 1) . Note how efficient gas drag partially erases the abrupt transition at 
(3 = 0.5 for the dynamical behaviour of small grains. 



measured bv lHeap et all <)2000|) from the HST/STIS scat- 
tered light images: about —4.5 for the northeast side of 
the disc and about —5.5 for the southwest side. It should 
be noted that a slope index close to -5 for the scattered 
profile beyon d the parent body region is in agreement with 
the results of iLecavelier et all l)l996|) who proved this to 
be a natural tendency for such systems. 

3.2. High gas density 

In this section we consider the high gas density case, with 
no = 1 6 cm~ 3 (case 1 in Tab. Ul, as inferred bv | Thi et alJ 
( 2001) and as required bv lBrandeker et all l)2004|) in order 
to have hydrogen to act as a "braking agent" for observed 
high-/? gaseous species. 

Differences with the gas free case remain negligible in 
the inner 100 AU region (FigQJ, where the flux is domi- 
nated by the biggest particles with (3 < 0.05 which are only 
moderately affected by gas drag (see FigEJl ■ Furthermore, 
it should be noted that the limiting value /Steady , sepa- 
rating the regime where gaseous friction leads to inward 
drift (/3 < /3 s teady) from the regime of outward drift (/? > 
/^steady), falls wi thin this size-range (see FigHJs). Using 
Equ. 27 and 28 of lTakeuchi fc Artvmowiczl (|200lh with the 
present gas density, one gets indeed /Steady = 0.034. 



In the r > 150 AU region, however, differences with the 
previous case are striking. The o btained brightness profile 
is clearly incompatible with the iHeap et al.l l|2000t) data, 



with F t 



400/100 



2.10 1 and p out 



-2.7. This is mainly 



due to the strong increase of the 0.05 < (3 < 0.5 parti- 
cles contribution. This population is strongly affected by 
gaseous friction and is progressively driven into the outer 
disc, hence flattening the radial profile, whereas in the gas 
free case only particles very close to the (3 — 0.5 limit have 
orbital eccentricities high enough to reach these regions 
(Fig. |2J). This outward drift has the additional effect of 
lengthening the collisional lifetime of these particles (be- 
cause of the radial dependency of the collision rate, cf. 
Equ.EJ). The contribution of the (3 > 0.5 particles is also 
increased with respect to the gas free case. For these par- 
ticles the effect of gas drag is on the contrary to slow down 
their radial escape, hence prolonging their stay in the disc 
(Fig.|21). Their contribution to the total flux remains nev- 
ertheless marginal. It is one order of magnitude smaller 
than the flux produced by the 0.05 < (3 < 0.5 particles 
between 100 and 400 AU (Fig.©. 
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Fig. 3. Simulated scattered light mid-plane profile in a 
gas free disc (solid line). The flux is arbitrary rescaled so 
that its maximum value equals f . The dotted line shows 
the contribution of the 0.005 < f3 < 0.05 population, the 
dashed line that of the 0.05 < j3 < 0.5 and the dash-dot 
line that of the 0.5 < < 25. The grey area shows the 
rang e of plausibl e profi les in the outer region, as deduced 
from lHeap et alJ <|2000j) observations for the SE and NW 
sides. It corresponds to the area bordered by r -4 5 and 



-5.5 



slopes derived from the flux value at 150 AU. 



Fig. 5. Global profile for different gas disc prescriptions. 
Solid line: gas free profile of Fig. [21 Then, from bottom 
to top of the profile value s at 500 AU: Solar abu ndance 
reconstructed gas disc, i.e. iBrandeker et "all l|2004h profile 
with no = 2.25 10 3 cm~ 3 , same profile with no = 10 4 cm~ 3 , 
IBrandeker et alJ ll2004fl profile with no = 10 6 cm -3 trun- 
cated at 150 AU. IBrandeker et all l)2004j) profile with n = 
10 5 cm -3 , Hayashi profile with a total disc mass of 40 M®. 
These models correspond respectively to cases 2, 3, 5, 4 
and 6 in Tabled 




30 40 50 60 70 80 90100 200 

radius (AU) 



300 400 500 



Fig. 4. Same as FigOU but for the high gas density disc 
(case 1 in Tabled . 



3.3. Exploring the gas density 

Taking the low gas den sity distribution reconstructed by 
IBrandeker et al. l|2004|) assuming solar abundances leads 
to a luminosity profile which is almost undistinguishable 
from the gas free one (Fig.[5J). Additional runs assuming 
the profile of Equ.|7|with varying no show that departures 
from the gas free cas e , and thus incompatibility with the 
observed iHeap et alJ l)2000(l luminosities, become signifi- 
cant for no > 10 4 cm~ 3 . We also performed an academic 



run with the high gas density distribution truncated be- 
yond 150 AU. But this leads to a luminosity profile much 
too high around 150-300 AU, mainly because of the accu- 
mulation of the 0.05 < (3 < 0.5 g rains in this region. This 
is an effect already observed bv iTakeuchi &: Artvmowic3 
i|200l|) for small grains in truncated gaseous discs. An ad- 
ditional run with a total gas mass equal to that of case 



2 (40 Mm) but assuming a solar type profile in r — 
( Hayash also leads to strong departures from the 

observed luminosity profile in the outer region. 

Note also that the dust luminosity profile in the inner 
100 AU region is almost independent from gas drag effects, 
at least for the density range explored here. 

3.4. Exploring the dust size distribution in the high gas 
density case 

All luminosity profiles sofar were computed assuming that 
dust particles follow the classical dA^ oc s~ 3 5 ds size dis- 
tribution, whic h is the prescription us ually taken in simi- 
lar studies (e.g. lAueereau et alJ . l200l|) . However, it should 
be stressed that "real" size distributions may significantly 
depart from this theoretical collisional equilibrium power 
law, especially for pa rticles close to the blow -out size (the 
reader might refer to lThebault et ail l)2003|) for a detailed 
discussion on this issue). 

Accurately estimating the size distribution exceeds the 
scope of the present work, but we have tested alternative 
size-distributions power laws in order to explore how the 
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Fig. 6. Profiles obtained for different dust size distribu- 
tions in the high gas density disc (case 1). dash-dot line: 
dN oc s~ 35 ds, dashed line: dN oc s _3 ds, dotted line: 



dN 



,-2.5 



ds. The solid line is the reference gas free pro- 



file of Fig.El 

final profile depends on this parameter. Figure|f)] clearly 
shows that reducing the power-law index of the size dis- 
tribution indeed leads to a drop of flux beyond 150 AU. 
However, the price to pay for this better fit in the outer 
disc is a significant and troublesome flux reduction in the 
100-150 AU region. This leads to profiles which are as bad 
a fit to observational data as in the dN oc s _3 5 ds case. 
This flux reduction could in principle be compensated by 
allowing the radial distribution of the parent bodies to 
extend further out. However, this would in turn lead to 
an increase of the flux beyond 150 AU, and yet departure 
from the good match originally obtained there, because 
the luminosity profile in these outer regions is mainly due 
to particles which originate from the highest parent body 
density region. 

4. Discussion 

The above results all tend to point towards an incom- 
patibility between high gas densities and dust luminosity 
profiles matching observations. One could wonder however 
if these discrepancies could not be eliminated by assuming 
pa rent bodies distr i bution s that differ from that proposed 
bv I Augereau et all 1 200lh in the context of a gas free sys- 
tem. In this case the problem would get hopelessly de- 
generated: for every gaseous disc considered there would 
always be an ad-hoc parent body distribution producing 
a good fit to the observations. 

We do not believe this to be the case because the sim- 
ulations show that parent body distributions cannot be 
arbitrarily chosen neither in the inner nor in the outer 
regions of the disc. In the inner (< 100 AU) region, the lu- 
minosity profile is dominated by big grains with f3 < 0.05 
(see Fig.0J). These particles are only weakly affected by gas 
drag and remain on nearly circular orbits matching those 



of the parent bodies producing them. It follows that, for 
the inner disc, a change of parent body distribution would 
automatically lead to departures from the best-fit lumi- 
nosity profile which cannot be compensated by gaseous 
friction effects, at least within the gas density range ex- 
plored in this paper (additional test runs have shown that 
significant departures from the gas free profile are found, 
in the inner disc, when no > 10 7 cm -3 , which corresponds 
to unrealistic large amounts of gas). 

Beyond 150 AU the problem gets more complex and 
there might in principle be an alternative parent body dis- 
tribution able to "bend" the outer disc profile towards the 
"right" one even for the high gas density case. However, it 
is easy to see that in order to achieve this, this alternative 
distribution should act as to decrease the outer disc lumi- 
nosity by more than an order of magnitude (Fig.^J). This 
cannot be done by adjusting the parent body distribution 
in the outer disc, since in the reference case this distri- 
bution is already truncated beyond 150 AU, so any alter- 
native distribution in the outer disc can only add objects 
there and result in an increase of the luminosity. Thus 
in order to lower the luminosity profile in the outer disc 
the only possibility is to change the parent body distribu- 
tion in the inner disc. But, as has just been shown, this 
cannot be done without in turn strongly modifying the 
luminosity profile in the inner regions (which would result 
in departing from a good fit of observations), thus solving 
one problem by creating another. A s a consequence, we 
do believe the lAueereau et alJ l)200l|) parent body profile 
to be a rather robust assumption. 

Alternatively, one may consider different dust size dis- 
tributions, for instance a shallower distribution in order to 
lower the flux in the outer regions where small grains dom- 
inate. However, Fig.EJshows that this doesn't work either. 
We cannot completely exclude the possibility that a fine- 
tuned fit might be found by adjusting the size distribu- 
tion at every distance from the star, but such a thorough 
parameter exploration is beyond the scope of the present 
work. However, even if such an ad-hoc fit should be found, 
it would have the disadvantage of being less generic. 

Another parameter space that was left unexplored 
here is the grain scattering properties. However, for 
dust radial distributions rela tively similar to the present 
ones, lAueereau et all l)200l|) has shown that assuming 
anisotropic scattering properties for the grains would not 
significantly affect the flux profile beyond ~80 AU, i. e. the 
region of interest for the present discussion. 

Under these assumptions we thus do believe that the 
present simulations rule out the presence of large amounts 
of gas around (3 Pictoris, espe cially in the outer parts o f the 
disc. Assuming the profile of iBrandeker et ahl (Hp04), the 
upper limit on the no value is ~ 10 4 cm~ 3 corresponding 
to a total gas mass of ~ 0.4M ffl . This corresponds to a 
gas to dust ratio smaller or approximatively equal to 1 at 
most. 
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5. Summary and conclusion 

We have numerically investigated the dynamics of dust 
particles suffering a gas drag force in the edge-on (3 Pictoris 
disc. From our simulations we were able to compare the- 
oretical scattered light profiles along the disc m i dplan e 
to the reference images obtained bv iHeat) et alJ l|2000|) . 
Several gas disc prescriptions and a broad range of grain 
sizes (from (3 = 0.005 to f3 = 25) have been explored. We 
assumed that dust particles are produced by parent bod- 
ies fol lowing the distribution derived by lAugereau et alJ 
l)200l[) . an assumptiom which our simulations prove to be 
relatively robust. 

Our main result is that high density gas discs always 
lead to dust distributions whose luminosity profiles are 
strongly incompatible with the observations in the outer 
regions of the disc (beyond 150 AU), while the surface 
brightness profile in the inner regions is poorly sensitive 
to the effect of g as drag on the grains. A ssuming the gas 
density profile of iBrandeker et ahl (|2004^) . gas discs with 
total mass above ~ 0.4M $ seem to be ruled out. This 
value is consistent wit h the H 2 upper limit obtained by 
iLecaveli ere^d] (|2f]0ll) and with the total gas mass esti- 
mated by lKampet alJ l|2003h . It is also consi s tent w ith the 
gas mass upper limit found bv IChen et al.1 l)2004[) based 
on recent Spitzer mid-infrared observations. Unless the 
gas disc is surprisingly very extended so that the dust 
dynamics is not affected by gas drag within the first hun- 
dreds of AU from the star, our results tend to rule out 
the much highe r gas mass estimate (50M®) derived by 
iThi et af] l|200l|) . However, this leaves open the puzzling 
problem of high f3 ions observed on bound orbits, since 
it also rules out the gaseous braking agent hypothesis for 
w hich at least 40Mm of h ydrogen are required according 
to IBrandeker et all 1 20041) . Our results suggest that other 
scenarios, like braking by other chemical species or by 
magnetic field, have definitely to be investigated. This may 
lead to a better estimate of the mass of gaseous species 
insensitive t o stellar radiation pres sure required to stop 
high /3 ions i Beust & Valiroil feoOS). The approach devel- 
opped in the present paper could also in the future be 
applied to similar debris disc systems, in particular the 
recently discove red disc around A UMic, which is coeval 
with /3 Pictoris l|Kalas et allr2004h . 
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